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ABSTRACT 

HAZX  is  a  software  tool  for  assessing  explosive  hazards  when  the  Quantity-Distance  (Q-D)  safe  separation 
requirements  are  violated.  This  paper  presents  the  development  of  DebrisHaz,  a  sub-module  in  HAZX  for 
performing  detailed,  site  specific  fragmentation  hazard  analyses.  Features  of  this  module  include:  a)  3D  modeling 
of  site  terrain  and  structures;  b)  random  sampling  of  fragments  from  a  user- specified  fragment  source;  c) 
computation  of  fragment  trajectories  with  drag  correction;  d)  physics-based  3D  simulation  of  fragment  impact, 
bounce  and  roll;  e)  use  of  trajectory-normal  method  for  calculating  fragment  density;  f)  capability  to  track  fragment 
impacts  with  a  3D  structure’s  wall  and  roof  surfaces  using  collision  detection  algorithms;  g)  calculation  of  the 
expected  number  of  penetrations  into  a  building.  The  trajectory  calculator  has  been  verified  with  a  range  of 
established  computer  programs,  and  the  bounce  model  has  been  compared  against  field  test  data  with  good 
agreement.  The  tool  has  been  used  to  perform  detailed,  site-specific  fragment  risk  analyses  that  account  for  unique 
fragmentation  characteristics,  uneven/non-uniform  terrain,  and  the  shielding  of  people  or  structures  by  natural 
geographical  features  or  by  other  structures. 


INTRODUCTION 

Quantity-Distance  (Q-D)  analysis  is  usually  the  first  step  in  the  process  of  gaining  explosive  site  plant  approval.  If 
a  site  plan  is  in  violation  of  any  of  the  Q-D  requirements,  additional  supplemental  hazard/risk  analyses  may  be 
performed.  These  supplemental  analyses  usually  take  into  account  more  detailed,  site  specific  information  that  is 
outside  the  consideration  in  a  Q-D  analysis.  Examples  of  such  site  specific  information  include  unique 
fragmentation  characteristics,  uneven/non-uniform  terrain,  and  the  shielding  of  people  by  natural  geographical 
features,  man-made  barriers  or  by  other  structures.  DebrisHaz,  one  of  the  modules  in  HAZX  (Chrostowski,  2010), 
is  a  tool  developed  to  perform  detailed,  site-specific  debris  hazard  analysis. 

In  a  DebrisHaz  analysis,  3 -dimensional  models  are  built  for  the  terrain  and  buildings  in  a  graphical  user  interface. 
Fragment  groups  are  developed  for  potential  explosion  scenarios.  Fragments  are  randomly  sampled  from  the 
fragment  groups;  their  drag-corrected  trajectories  are  tracked  and  their  impacts  with  the  terrain  and  buildings  are 
determined.  In  the  case  of  a  terrain  impact,  the  translational  and  rotational  bounce  velocities  are  determined  with  a 
physics-based  bounce  analysis,  and  the  trajectory  tracking  continues  with  the  new  velocities.  In  the  case  of  a 
building  impact,  the  potential  of  wall  or  roof  penetration  is  evaluated.  If  the  trajectory  of  a  fragment  is  below  the 
height  of  a  typical  person,  its  contribution  to  the  pseudo  trajectory  normal  (PTN)  fragment  density  is  calculated  and 
tallied. 
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Outputs  from  the  program  are  the  PTN  fragment  density  distribution  over  the  terrain  and  the  expected  number  of 
penetrations  for  each  building.  The  PTN  fragment  density  can  be  contoured  to  determine  the  hazardous  fragment 
density  area  and  distance.  The  building  penetration  data  can  be  used  for  further  risk  analysis. 

FRAGMENT  GROUP 

A  fragment  group  consists  of  one  or  more  fragments  with  statistically  identical  characteristics,  which  include  the 
mass,  size,  shape,  aerodynamic  properties,  initial  location  and  the  explosion-induced  velocity  of  the  fragments. 
Each  fragment  group  is  defined  in  a  reference  frame  that  is  usually  different  from  the  global  coordinate  system  in 
which  the  buildings  and  terrain  are  defined.  The  reference  frame  for  a  fragment  group  can  be  stationary  or  moving. 
The  latter  allows,  for  example,  the  program  to  simulate  explosion  scenarios  due  to  a  space  launch  failure.  Each 
simulation  can  contain  one  or  more  fragment  groups. 

During  a  Monte  Carlo  (MC)  simulation,  fragments  are  randomly  sampled  from  each  fragment  group.  All  vector 
parameters  are  converted  from  the  local  reference  frame  to  the  global  coordinate  system,  and  the  velocity  of  a 
moving  reference  frame  is  added  to  the  explosion-induced  velocities  of  individual  fragments. 

FRAGMENT  BALLISTIC  TRAJECTORY 

The  translation  and  rotation  of  a  ballistic  fragment  are  described  by  the  Newton-Euler  equations: 

m*c(t)  =  'Ef  (!) 

Icb(y)  +  to(7)xIco(y)  =  (2) 

where  m  is  the  fragment  mass,  I  is  the  fragment  moment  of  inertia  matrix,  r,  is  the  translational  position  vector  of 
the  fragment  center  of  mass,  CO  is  the  fragment  angular  velocity,  is  the  sum  of  all  forces  acting  on  the 

fragment,  and  ^  T  is  the  sum  of  all  moments  acting  on  the  fragment.  The  forces  acting  on  the  fragment  include 

the  gravity  and  the  drag.  In  a  global  coordinate  system  with  the  z-axis  pointed  upward,  the  vector  of  the  gravity  can 
be  written  as 


The  drag  can  be  written  as 
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fd  =  ~\cdAPair  ||v||  v  =  ~^CdAPair  ||rc||rc .  (4) 

where  CD  is  the  coefficient  of  drag  and  A  is  the  presented  area  of  the  fragment.  Substituting  Equations  (3)  and  (4) 
into  (1)  gives 
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(5) 


The  fragment  is  assumed  to  be  free  of  any  moments  while  it  is  in  free  flight,  and  so  the  right-hand  side  term  in 
Equation  (2)  is  zero.  The  initial  translational  velocity  of  a  fragment  is  the  explosion-induced  velocity  plus  the 


velocity  of  the  moving  frame.  The  initial  rotational  velocity  of  a  fragment  is  assumed  to  be  zero,  although  the  first 
impact  with  the  terrain  will  usually  introduce  rotation.  The  equations  for  translation  and  rotation  are  integrated 
numerically. 


FRAGMENT  BOUNCE  ANALYSIS 


Bounce  Equations 


The  bounce  analysis  determines  the  translational  and  rotational  velocities  of  a  fragment  after  it  strikes  the  ground 
surface.  The  analysis  is  based  on  an  impulse  collision  response  algorithm  described  in  (Mirtich,  1996).  The  basic 
assumptions  of  the  analysis  are:  (1)  the  deformation  of  the  fragment  and  the  impacted  surface  are  negligibly  small 
so  that  they  can  be  considered  rigid;  (2)  Stronge’s  hypothesis  applies;  and  (3)  Coulomb’s  friction  law  holds. 

The  first  assumption  implies  that  the  duration  of  the  impact  is  infinitesimal  and  the  force  on  the  fragment  is 
impulsive  and  instantaneously  changes  the  velocity  of  the  fragment.  Since  the  fragment  velocities  are  finite,  the 
position  of  the  fragment  remains  constant  during  the  impact.  Forces  with  finite  magnitudes  (such  as  gravity)  have 
no  effect  over  infinitesimal  intervals,  and  so  can  be  ignored  in  the  bounce  analysis. 

While  the  first  assumption  idealizes  the  fragment  and  the  impacted  surface,  in  reality,  there  is  always  deformation 
associated  with  a  force.  The  entire  impact  can  be  divided  into  two  phases:  a  compression  phase  and  a  restitution 
phase.  During  the  compression  phase,  the  fragment  pushes  into  the  surface  until  the  normal  component  of  the 
fragment  velocity  at  the  contact  point  becomes  zero;  during  the  restitution  phase,  the  fragment  is  pushed  out  by  the 

restitution  force.  Let  Wz  ( tc )  and  Wz  (Vr)  be  the  work  done  by  the  normal  component  of  impacted  surface’s 

resistance  during  the  compression  phase  and  the  restitution  phase,  respectively,  then  Stronge’s  hypothesis  states  that 

W,(tr)  =  -e‘Wz(tc)  (6) 

where  e  is  the  coefficient  of  restitution,  and  the  negative  sign  reflects  the  fact  that  the  work  done  during  the  two 
phases  have  opposite  signs. 

In  simple  impact  situations,  such  as  for  a  spherical  fragment,  Stronge’s  hypothesis  reduces  to  Newton’s  impact  law: 


The  Coulomb  friction  law  can  be  expressed  as 

ft  =-#/>, /hi  for  u,*°  (8) 

\%\\^Mfn  for u,  =  0 

where  fl  is  the  coefficient  of  friction,  Ut  is  the  tangential  component  of  the  fragment  velocity  at  the  impact  point, 
fn  is  the  magnitude  of  the  normal  compression  force,  and  ft  is  the  tangential  friction  force.  With  the  non- 
impulsive  forces  dropped,  the  Newton-Euler  equations  for  the  fragment  in  an  impact  can  be  written  as 


(9) 


f  (7)  =  ma(?) 

r  xf  (V)  =  Id)(V)  +  ®(V)x(i(o(V)) 


(10) 


where  m  and  I  are  the  fragment  mass  and  moment  of  inertia  matrix,  a(7)  is  the  translational  acceleration  at  the  center 
of  mass,  (0  (7)  is  the  angular  velocity,  f(f)  is  the  impact  force,  and  r  is  the  offset  vector  from  the  center  of  mass  to 

the  impact  point,  as  shown  in  Figure  1.  It  is  noted  that  the  ground  surface  needs  not  to  be  horizontal  and  the  above 
equations  are  given  in  the  body  frame,  a  reference  frame  attached  to  the  fragment. 


Figure  1.  Fragment  impact  on  a  surface 


The  impact  force  and  the  resulting  translational  and  rotational  accelerations  dominate  the  above  equations,  so  the 
last  term  of  Equation  (10)  can  be  dropped.  With  these  simplifications,  the  above  equations  can  be  formally 
integrated  over  time  to  give 


p(7)  =  mAv(7)  (11) 

rxp(f)  =  IAco(f)  (12) 

where  Av(^)  and  Aco(V)  are  the  changes  in  the  translational  and  angular  velocities,  and  p(f)  is  the  impact 
impulse: 


t 

p(/)  =  Jf(r)/r.  (13) 

0 

Time  t  has  been  used  in  the  above  formal  derivation.  For  rigid  body  impact,  the  impact  duration  approaches  zero 
and  the  impact  forces  become  infinite.  To  avoid  difficulties  with  these  conditions,  the  impulse  and  velocities  can  be 
expressed  as  functions  of  a  new  impact  parameter  that  monotonically  increases  during  the  course  of  the  impact. 

Denoting  the  new  impact  parameter  by  X  ,  Equations  (11)  and  (12)  can  be  written  as 

AvM  =  ~pW  (14) 

m 

Aco(/l)  =  I-1  [r  . 


The  contact  point  velocity  u  is  given  by 


(15) 


u(/l)  =  v(/l)  +  o>(/l)xr 


(16) 


and  the  change  of  u  can  be  written  as 


Au(/l)  =  Av(/l)  +  A«(/l)xr .  (17) 

In  the  following  derivation,  the  impact  parameter  is  dropped  from  the  variables  unless  ambiguity  arises. 
Substituting  Equations  (14)  and  (15)  into  (17)  gives 
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and  ru  r2,  and  r3  are  the  components  of  vector  r.  Let 

K  =  — l-rlr 

m 


then  Equation  (18)  becomes 


Au  =  Kp . 


(18) 


(19) 


(20) 


(21) 


Matrix  K  is  constant,  symmetric  and  positive  definite  (Mirtich,  1996).  The  differential  form  of  Equation  (21)  can 
be  written  as 


d\X  _  ^  t/p 

d/ l  d/ l 


(22) 


The  above  derivation  is  based  on  the  Newton-Euler  equations  established  in  the  body  frame.  It  is  more  convenient 
to  carry  out  the  bounce  analysis  on  a  reference  frame  whose  z-axis  is  aligned  with  the  normal  of  the  impacted 
surface  and  whose  x-  and  y-axis  fall  on  the  impact  surface.  This  reference  frame  is  denoted  the  space  frame. 


Let  R  be  the  transformation  from  the  body  frame  to  the  space  frame.  Then  in  the  space  frame  the  vectors  u  and  p 
can  be  written  as 
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ii 
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(23) 

II 

(24) 

Solving  for  u  and  p  from  Equations  (23)  and  (24)  and  substituting  them  into  Equation  (22)  gives 


where 


K  =  RKRr . 

Since  K  is  positive  definite,  K  is  also  positive  definite,  and  Equation  (25)  can  also  be  written  as 

t/p  _j£-i  ^ 

d/ l  d/ l 

In  the  case  where  the  tangential  velocity  is  non-zero,  that  is, 

yl%+%  >0 


(26) 


(27) 


(28) 


the  Coulomb  friction  law  gives  the  force  at  the  contact  point  as 


where 
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(30) 


If  ^  u l  +U2y  —  0  ,  then  ^  is  undefined.  This  situation  is  termed  sticking  and  requires  special  treatment.  Sticking  is 

divided  into  stable  sticking  and  unstable  sticking.  In  stable  sticking,  the  tangential  velocity  becomes  zero  and 
remains  zero  from  thereon.  In  non-stable  sticking,  the  tangential  velocity  becomes  zero  only  momentarily  and 
sliding  resumes  after  that. 

The  integration  of  the  bounce  equations  is  divided  into  the  compression  phase  and  the  restitution  phase.  In  the 
compression  phase  integration,  the  normal  component  of  the  fragment  velocity  Uz  at  the  contact  point  is  used  as  the 

independent  variable,  and  the  tangent  components  of  the  contact  point  velocity,  Ux  and  Uy  ,  and  the  work  Wz  done 

by  the  normal  component  of  the  resistance  force  are  the  dependent  variables.  The  integration  is  carried  out  until  uz 
becomes  zero. 


In  the  restitution  phase  integration,  Wz  is  used  as  the  independent  variable  and  the  fragment  velocity  at  the  contact 
point  is  the  dependent  variable.  The  integration  is  carried  out  until  Wz  reaches  a  value  determined  from  Stronge’s 
hypothesis.  The  reader  is  referred  to  (Mirtich,  1996)  for  details  of  the  integration  process. 

The  outcome  of  the  restitution  phase  integration  is  the  bounce  velocity  U  at  the  contact  point.  The  velocity 
change,  Au  ,  is  the  difference  between  the  bounce  velocity  and  the  initial  velocity.  The  total  impulse  that  acts  on 
the  fragment  during  the  impact  can  be  calculated  from  Equation  (27)  as 

p  =  K  ‘Au 

in  the  space  frame,  or 

p  =  Rrp  =  RrK ‘Au.  (32) 


(31) 


in  the  body  frame. 

The  change  in  the  angular  velocity  is  then  calculated  from  (15)  as 

Aco  =  I“‘(rxp).  (33) 

The  change  in  the  center  of  mass  velocity  can  be  calculated  from  Equation  (14)  as 

Av  =  Tp  (34) 

m 

in  the  body  frame.  With  Ao  and  Av  known,  the  exit  angular  velocity  and  translational  velocity  can  be  calculated: 


(tf )  =  to  (f0 )  +  Ato 

(35) 

II 
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£ 

(36) 

The  exit  angular  velocity  and  translational  velocity  are  used  as  the  initial  conditions  for  the  next  phase  of  airborne 
trajectory  of  the  fragment. 

Comparison  with  Experimental  Data 

Comparison  of  the  bounce  model  was  made  with  the  experimental  data  presented  in  (Fletcher  &  Bowen,  1968) 
based  on  the  ending  bounce  distance  of  concrete  blocks.  In  the  experiment  more  than  500  objects  were  dropped  on 
smooth  flat  ground  (a  graded  airstrip  on  an  alluvial  plain)  from  a  truck  traveling  at  speeds  ranging  from  10  to  60 
mph  and  the  distances  traveled  by  the  objects  were  measured.  The  objects  included  stones  (1.3 — 116  lb),  ordinary 
concrete  building  blocks,  concrete  filled  building  blocks,  and  several  types  of  animals.  Our  comparison  study  was 
made  only  with  the  filled  concrete  building  blocks  since  their  weights  (55.0-56.5  lb),  shape  (cuboid),  and  size 
(7.5x7.7x15.5  inches)  were  well  defined. 

In  the  analytic  simulation,  the  blocks  were  dropped  in  random  orientations  at  the  height  of  the  truck  bed  (2.6  feet) 
and  with  a  horizontal  velocity  ranging  from  15  ft/s  to  100  ft/s,  mimicking  the  test  conditions.  The  coefficient  of 
restitution  and  the  coefficient  of  friction  were  not  available  from  the  test.  Instead,  the  restitution  and  friction 
coefficients  used  in  ACTA’s  range  safety  software  were  employed  in  the  simulation.  The  coefficient  of  restitution 
is  a  function  of  the  normal  impact  speed,  as  shown  in  Figure  2.  The  coefficients  of  friction  are  constant,  with  0.9 
for  soft  surface  and  0.5  for  hard  surface. 


Figure  2.  Coefficient  of  Restitution. 

Figure  3  shows  sample  fragment  trajectories.  All  the  trajectories  in  the  figures  correspond  to  the  same  initial 
fragment  height  and  velocity.  The  initial  orientation  of  the  fragment  is  randomly  sampled  and  is  responsible  for  the 
variations  in  the  trajectories  after  the  first  impact.  It  is  noted  that  many  of  the  fragments  reach  a  height  above  the 
initial  drop  height  of  2.6  feet.  This  is  due  to  the  conversion  of  the  momentum  associated  with  the  horizontal 
velocity  into  vertical  momentum,  and  is  one  of  the  differences  between  a  spherical  object  and  a  non-spherical 
object. 


Figure  3.  Sample  fragment  trajectories  (height  vs.  downrange  &  cross  range  vs.  down 

range). 


Error!  Reference  source  not  found,  compares  the  test  and  simulated  fragment  travel  distances.  In  the  figure,  the 
solid  and  hollow  squares  represent  the  mean  test  distances  for  the  solid  and  the  hollow  blocks  respectively.  The 
regression  line  and  the  lines  for  the  ±1  standard  error  for  the  test  data  are  plotted  in  solid  and  dashed  lines.  The 
simulation  results  with  soft  and  hard  surfaces  are  plotted  in  error  bars  (mean  ±1  standard  error).  As  one  can  see,  the 
predicted  distances  are  higher  than  the  test  results  in  most  cases  but  are  still  within  the  range  of  error  of  the  test. 


The  over  prediction  with  the  hard  surface  is  more  pronounced  than  with  the  soft  surface,  which  is  expected 
considering  the  fact  that  the  hard  surface  model  is  intended  for  concrete  pavement.  Figure  5  shows  a  sequence  of 
3D  snapshots  from  a  bounce  simulation 


Figure  4.  Comparison  with  experimental  stopping  distance. 

PSEUDO  TRAJECTORY  NORMAL  FRAGMENT  DENSITY 

Fragment  density  (number  of  fragment  impacts  per  unit  area)  is  a  common  used  quantity  to  measure  the  hazard  of 
fragment  impact.  Depending  on  the  situations,  the  fragment  density  can  be  calculated  with  a  number  of  different 
methods,  and  the  results  from  these  methods  may  not  be  comparable.  In  field  tests,  for  example,  fragments  in 
gridded  areas  are  picked  up  after  they  have  stopped.  In  simple  numerical  simulations,  the  number  of  fragment 
impacts  is  usually  counted.  Neither  method  is  able  to  distinguish  the  difference  between  a  high  angle  fragment  that 
falls  almost  vertically  and  a  low  angle  fragment  that  has  a  flat,  near  ground  trajectory.  In  both  of  the  two  methods,  a 
fragment  flying  past  a  grid  area  3  feet  above  ground  will  have  no  contribution  to  the  fragment  density  for  the  grid 
area,  even  though  a  person  standing  inside  the  grid  area  can  potentially  be  hit  by  the  fragment.  The  pseudo 
trajectory  normal  (PTN)  fragment  density  is  an  alternative  method  that  is  used  to  correct  these  problems. 

The  PTN  fragment  density  accounts  for  the  height  Hr  of  the  impact  receiver  of  interest,  and  its  calculation  is  based 
on  the  intersection  of  the  fragment  trajectory  with  an  associated  volume  generated  by  vertically  extruding  a  grid 
area  on  the  ground  surface  to  a  height  equal  to  Hn  as  shown  in  Error!  Reference  source  not  found,  for  a  triangular 
grid  area.  The  PTN  fragment  density  for  a  particular  ground  triangle  is  defined  as 

Dptn  =  ^ 

^ PTN 


(37) 
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Figure  5.  A  sequence  of  frames  from  a  fragment  bounce  simulation. 


Figure  6.  Pseudo  Trajectory  Normal  Fragment  Density  Calculation. 


where  the  summation  is  over  all  fragments  whose  trajectories  intersect  with  the  associated  volume  of  the  ground 
triangle,  and  Apm  is  the  PTN-projected  area  of  the  associated  volume  with  respect  to  a  fragment  trajectory.  The 
PTN-projected  area  for  an  associated  volume  is  the  sum  of  the  PTN-projected  areas  of  the  faces  of  the  prism: 

(38) 

where  the  summation  is  over  all  faces  of  the  prism,  and 
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'4ni  •  v 
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if  n,  •  x f  <  0 
if  n;  •  \f  >  0 


(39) 


where  At  and  nz  are  the  area  and  unit  (outward)  normal  vector  of  z-th  face  of  the  prism,  respectively,  and  is  the 
fragment  velocity  vector. 

It  is  noted  that,  if  all  the  fragments  fall  vertically,  the  PTN  method  gives  the  same  fragment  density  as  dividing  the 
number  of  fragments  that  hit  a  ground  triangle  by  the  area  of  the  triangle.  The  two  methods  diverge  for  fragments 
falling  in  slanted  trajectories.  The  more  the  trajectories  are  slanted,  the  larger  the  difference.  As  the  trajectories 
become  horizontal,  the  number  of  fragment  impacts  on  the  ground  triangle  becomes  zero,  but  the  PTN  still  gives  a 
reasonable  fragment  density  that  can  be  used  to  evaluate  the  hazard  of  fragment  impact  on  exposed  people. 

Figure  7  shows  the  PTN  fragment  density  due  to  a  simulated  explosion  of  a  jet  assisted  takeoff  solid  rocket  motor 
used  in  the  launch  of  target  planes.  The  motor  is  tilted  up  and  pointed  to  the  east  to  northeast  direction.  The 
preferential  dispersion  of  the  fragments  in  the  radial,  axial  forward  and  axial  rearward  directions  give  the  butterfly- 
shaped  fragment  density  contour. 
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Figure  7.  PTN  fragment  density  due  to  a  simulated  explosion  of  a  jet  assisted  takeoff 

motor. 


FRAGMENT  IMPACT  ON  3D  OBJECTS 


The  PTN  fragment  density  discussed  above  is  appropriate  for  evaluating  the  fragment  impact  hazard  for  people 
since  the  height  of  an  average  person  is  used  as  the  characteristic  height  in  the  calculation.  While  the  same 
calculation  can  be  done  for  3 -dimensional  objects,  such  as  buildings  and  storage  tanks,  difficulties  arise  due  to  the 
potentially  wide  variation  in  the  heights  of  these  entities.  Additionally,  we  may  also  be  interested  in  the  penetration 
of  fragments  into  buildings,  which  is  affected  by  construction  of  the  walls  and  roof  of  the  buildings  as  well  as  the 
impact  angle.  Under  such  a  circumstance,  a  more  appropriate  approach  is  to  determine  the  impacts  on  each  3D 
object  directly  during  a  simulation.  To  do  this,  a  3D  object  is  modeled  with  triangle  mesh,  as  shown  in  Figure  8, 
and  the  impact  of  a  fragment  with  the  object  is  determined  through  collision  detection. 

Figure  9  shows  the  results  of  a  simulated  early  time  launch  failure.  The  trajectory  plot  includes  only  a  randomly 
selected  subset  of  all  the  trajectories  generated  in  the  simulation.  The  buildings  are  colored  by  the  expected  number 
of  fragment  penetrations,  with  blue  representing  no  penetration  and  red  the  highest  number  of  penetrations. 
Fragment  bounce  is  not  included  in  this  particular  simulation.  Figure  10  shows  a  simulation  of  the  1997  Delta  II 
explosion. 


Figure  9.  Fragment  impact  in  a  simulated  early  time  launch  failure. 


Figure  10.  Simulation  of  the  1997  Delta  II  explosion. 

SUMMARY 

DebrisHaz  performs  detailed  3D  fragment  throw  simulation,  computes  pseudo  trajectory  fragment  density  on 
meshed  terrain,  and  outputs  expected  numbers  of  fragment  impacts  and  penetrations  for  buildings.  It  provides  a 
useful  means  to  perform  site-specific  analysis  in  supplement  to  the  regular  Q-D  analysis,  taking  into  account  of  the 
information  unique  to  the  project.  The  program  is  still  under  active  development,  and  more  features  will  be  added 
to  it.  DebrisHaz  will  be  incorporated  as  a  module  in  the  HAZX  Hazard  Tool  to  simplify  its  use  by  end  users 


REFERENCES 


Chrostowski,  J.  D.,  et  al  (2010),  HAZX  Parti:  An  Explosion  Hazard  Assessment  Tool,  2010  Department  of  Defense 
Explosives  Safety  Board  Seminar,  Portland,  Oregon. 

Fletcher,  E.  R.,  &  Bowen,  I.  G.  (1968).  Blast-Induced  Translational  Effects.  Annals  of  NY  Academy  of  Sciences  , 
752  (1),  378-403. 

Mirtich,  B.  V.  (1996).  Impulse-Based  Dynamic  Simulation  of  Rigid  Body  Systems.  PhD  Thesis,  University  of 
California,  Berkeley,  CA. 


DebrisHaz 

3D  Fragment  Throw  Simulation 
to  Determine  Fragment  Fragment  Hazards 

and  Consequences 

Author/Co-Author 

Wenshui  Gan/Jon  Chrostowski 

ACTA,  Inc. 


2010  DDESB  Explosives  Safety  Seminar 

Portland,  OR 

July  15,  2010 


Radiaf  Fragment 
HFDD  =130  feet 


operations. 


principles  apply  to  every  possible 


DDESB. 


rices,  or  ACTA.  Inc.  cannot 
ponsible  for  any  mishap  « 
resulting  from  decisions 
d  by  this  software.  The  u 
inalysis  software.  HAZX  i 


laoOens^Bfl 


HAZX 


DDESB  2010  Seminar  -  DebrisHAZ 

DebrisHaz 


Performs  detailed,  site-specific  3D 
fragment  throw  simulations  to 
determine  their  hazards  and 
consequences  to  people  in  the  open 
and  to  building  occupants. 

DebrisHAZ  is  being  incorporated 
underneath  the  HAZX  Hazard  Tool 
(HRT) 

Features: 

-  Uses  3D  terrain  acquired  from 
GoogleEarth,  NIMA  DIED  national 
database  or  other  source(s). 

-  Extrudes  buildings  from  footprints  to 
their  respective  heights 

-  Assigns  construction  types  to  building 
walls,  roof,  windows. 

-  Performs  physics-based  trajectory  and 
fragment  bounce  analyses. 

-  Computes  PTN  fragment  density. 

-  Computes  expected  numbers  of 
fragment  impacts  &  penetrations  of 
buildings. 
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Fragments  and  Fragment  Groups 


•  DebrisHaz  considers  3D  fragments 
and  point  mass  fragments. 

-  Shapes  include  spheres,  cuboids, 
ellipsoids,  rods 

•  A  fragment  group  contains 
fragments  with  statistically  identical 
properties  (random  but  w/  same  mean 
and  distribution  on  size,  mass,  velocity, 
etc.) 

•  Each  fragment  group  is  defined  in  a 
local  reference  frame,  which  can  be 
stationary  or  moving  (as  in  the 
explosion  of  an  in-flight  rocket). 

•  Fragments  are  randomly  sampled 
from  fragment  groups  and  their 
interactions  with  the  terrain  and 
buildings  are  tracked. 
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Free-Flight  Trajectory 


Newton-Euler  equations 


mvc  (,)  =  £f 

Id)(r)+co(^)x  to(f)=yT 


The  integration  of  the  N-E  equations  is  carried  out 
numerically  until  the  fragment  impacts  the  ground, 
at  which  point  the  calculation  is  switched  to  a  bounce 
analysis. 
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Fragment  Bounce  Analysis 


*  Basic  Assumptions 

-  A  fragment  has  shape  and  size 

-  Deformation  of  fragment  and  impacted  surface  is 
negligible. 

-  Stronge’s  restitution  hypothesis  applies 

•  Work  done  by  restitution  force  is  a  fraction  of  the  work 
done  by  the  compression  force. 

•  For  spherical  fragment,  this  reduces  to  Newton’s  impact 
law. 

-  Coulomb’s  friction  law  applies. 
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Fragment  Bounce  Analysis 


•  Mass 

•  Moment  of  inertia 

•  Forces  during  impact: 

-  Normal  reaction,  N 

-  Friction,  F 

•  Velocities 

-  Translational  velocity  at 
C.G. 

-  Rotational  velocity 
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Fragment  Bounce  Analysis 


Newton-Euler  equations: 
f  (?)  =  mv(?) 

rxf  (?)  =  Icb(?)+c)(?)x(lco(?)) 


Contact  point  velocity: 
u  =  v+coxr 

Contact  force: 


f  = 


u2  +u2 

x  y 


for  u2  +  ul  >  0 

x 


y 
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Fragment  Bounce  Analysis 


After  applying  the  assumption  of  impulsive  impact: 

p  =  mAv 

r  x  p  =  IAco 

Au  =  Ay  +  Aco  x  r 
p  =  fAt 

Integration  of  the  above  equations  give  the  bounce 
velocities. 
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Fragment  Bounce  Test  Comparison 


•  Test  data  by  Fletcher  &  Bowen  (1968), 

-  7.5”x7.5”x15.5”  hollow  &  solid  concrete  blocks 

-  Dropped  from  the  bed  of  a  truck  moving  at  speeds 
ranging  from  10  to  60  mph  on  a  graded  airstrip 

-  Travel  distance  measured 

•  Simulation 

-  Solid  concrete  block  modeled  as  cuboid 

-  Drop  height  2.6  ft  (height  of  truck  bed) 

-  Horizontal  velocity  15  ft/s  to  100  ft/s 


Vf 
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Fragment  Bounce  Test  Comparison 


Restitution 
coefficients  were  not 
available  from  test. 
Instead,  data  in 
ACTA’s  range  safety 
software  were  used. 

Coefficient  of  friction 
is  0.9  for  soft  surface 
and  0.5  for  hard 
surface. 

Simulations  were 
done  for  both  hard 
and  soft  surfaces. 
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Fragment  Bounce  Simulation  Samples 


Sample  fragment  trajectories  with  bounce 


Side  view 


Plan  view 


Note:  All  fragments  in  the  plots  have  the  same  initial 
velocity  and  height  but  their  takeoff  orientations  are  random. 
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Snapshots  of  the  3D  Bounce  Model 
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Stopping  Distance  (ft) 
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Fragment  Bounce  Test/Simulation  Comparison 
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Pseudo  Trajectory  Normal  Fragment  Density 


•  Hazardous  fragment  density  is 
used  in  determining  equivalent 
protection. 

•  Different  methods  exist  for  the 
calculation  of  fragment  density. 

-  Some  methods  may 
underestimate  the  real 
impact  risk  for  low-flying 
fragments  with  flat 
trajectories. 

-  Calculation  based  on  PTN 
corrects  this  problem. 


AwmV 
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Pseudo  Trajectory  Normal  Fragment  Density 


In  DebrisHaz  the  ground 
surface  is  meshed  with 
triangles. 

The  triangles  are 
extruded  to  the  height  of 
a  typical  person. 


Each  hazardous 
fragment  whose 
trajectory  intersects  with 
the  extruded  volume  of  a 
triangle  contributes  to 
the  fragment  density  of 
the  triangle. 
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Pseudo  Trajectory  Normal  Fragment  Density 


The  contribution  of  a 
ragment  to  the  PTN 
ragment  density  is 
the  inverse  of  the 
area  resulting  from 
projecting  the 
extruded  vo  ume 
onto  a  plane  that  is 
normal  to  the 
trajectory  of  the 
fragment. 
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Example  of  PTN  Fragment  Density  Calculations 


Explosives  Std 
HFDD  =  200  feet 


Radial  Fragment 
HFDD  =130  feet 


Simulated  on-rail  explosion  of  a  winged  target  vehicle  with  a 
solid  rocket  motor  assist  at  takeoff. 

Asymmetric  shape  of  the  hazardous  fragment  density  contour 
results  from  preferential  fragment  flight  out  pattern. 
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Fragment  Impact  on  Buildings 


•  PTN  fragment  density  is  suitable  for 
evaluating  the  risk  of  fragment  impact 
on  people  in  the  open. 

-  Intervening  buildings  and  other  objects 
may  provide  protection 

•  It  is  not  suitable  for  predicting  building 
impacts  and  consequences  to 
occupants: 

-  Wide  variation  in  building  heights 

-  Wide  variation  in  impact  resistance,  both 
within  a  building  (wall/roof)  and  among 
different  buildings. 

-  Blocking  of  fragments  by  one  building  on 
another 


•  The  DebrisHAZ  physics-based  tool 
addresses  all  these  issues 
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Fragment  Impact  on  Buildings 
(This  process  is  being  automated  in  HAZX) 


3D  models  are  created 
using  Google  Sketchup™. 

-  Footprint  are  extruded 
upward  to  height 

-  Surfaces  are  meshed 

Structural  properties  are 
assigned  to  walls,  roof  & 
windows. 

-  Properties  used  to 
determine  kinetic 
energy  of  penetration 

DebrisHaz  determines 
fragment  impacts  on 
building  surfaces  &  uses 
KE  &  frag  area  to  predict 
whether  it  penetrates  or  not 


*  Residual  fragment  KE  is  used 
to  predict  injury  to  building 
occupants 
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Examples  of  DebrisHaz  Analyses 
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Examples  of  DebrisHaz  Analyses 


1997  T+14  explosion 
of  a  Delta  II 
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Summary 


•  DebrisHaz  is  used  to  perform  detailed  3D  fragment  throw 
simulations  when  site-specific  siting  issues  are  important 

-  Tracks  trajectories  of  all  fragments  &  calculates  intersects 
with  other  objects  (e.g.,  bldgs,  barriers)  and  the  ground 
surface 

-  Considers  terrain  and  the  bounce  of  fragments 

-  Performs  a  PTN  analysis  to  determine  hazardous  fragment 
density 

-  Calculates  impacts  and  penetrations  of  buildings  and  effects 

_ on  occupants _ 

•  DebrisHaz  is  being  incorporated  underneath  HAZX  in 
order  to  automate  the  entire  process 

-  Users  will  be  able  to  perform  detailed  fragment  throw 
analyses  on  their  own 


